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The by now well-established scalar-isoscalar resonance /o(500) (the a meson) seems potentially 
relevant in the evaluation of thermodynamic quantities of a hadronic gas, since its mass is low. 
However, we recall that its contribution to isospin-averaged observables is, to a surprising accuracy, 
canceled by the repulsion from the pion-pion scalar-isotensor channel. As a result, in practice one 
should not incorporate /o(500) in standard hadronic resonance-gas models for studies of isospin 
averaged quantities. In our analysis we use the formalism of the virial expansion, which allows 
one to calculate the thermal properties of an interacting hadron gas in terms of derivatives of the 
scattering phase shifts, hence in a model-independent way directly from experimentally accessible 
quantities. A similar cancellation mechanism occurs for the scalar kaonic interactions between the 
7 = 1/2 channel (containing the alleged 77o(800) or the k meson) and the I = 3/2 channel. 
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I. INTRODUCTION 

Nowadays, there is no doubt about the existence of the 
scalar-isoscalar resonance /o(500) (traditionally known 
as the a meson). This very wide resonance is best de¬ 
scribed through its pole position, which according to the 
Particle Data Group (PDG) lies in a rather conservative 
range of M a -iT a /2 = (400 - 550)-z(200-350) []]. Such 
a low range for the mass was first reported in the 2012 
version of the PDG tables [3], while previous editions 
quoted a very broad range 400-1200 MeV. Detailed stud¬ 
ies of this resonance led to much smaller uncertainties: 
in Ref. [3] the pole at (400 ± 6 ^ 3 g) — z(278 ± 6 + 43 ) is ob¬ 
tained, while in Ref. [4HZ] the value (4571(g) — z(279l( 1 ) 
is quoted (for a compilation of all estimates and numer¬ 
ous references, see the mini-review ‘Note on scalar mesons 
below 2 GeV’ in Ref. HD- 

The question that we clarify in this paper is the role of 
the resonance /o(500) in thermal models for relativistic 
heavy-ion collisions (for reviews see, e.g., Refs. pH [9] and 
references therein). Naively, one is tempted to include it 
as a usual but broad Breit-Wigner resonance, as was done 
for instance in Refs. nam. Notice that the pole mass 
of /o(500) is light (even lighter than the kaon), hence it 
might have a non-negligible effect on observables in more 
accurate studies. Indeed, in Ref. El the effect of /o(500) 
was found to increase the pion yield by 3.5%. 

However, broad resonances should be treated with 
care. Moreover, there are also repulsive channels in the 
hadronic gas which must be taken into account on equal 
footing. The appropriate framework here is the virial 
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expansion (see, e.g., Ref. [13]). There, one uses only sta¬ 
ble particles (pions, kaons, nucleons, ...) as degrees of 
freedom, and their interactions are systematically incor¬ 
porated in thermodynamics via coefficients of the virial 
expansion. The second term of this expansion, corre¬ 
sponding to the 2 —> 2 reactions, is straightforward to 
include, as it involves the phase shifts accurately known 
from experiment. Within this approach, the attraction 
due to the pole of /o(500) is encoded in the pion-pion 
phase shift in the isoscalar-scalar (I = J = 0) channel. 
In principle, this attraction could generate non-negligible 
effects to thermodynamic properties. However, there is 
also repulsion from the isotensor-scalar (I = 0, J = 2) 
channel. It turns out that this repulsion, whose effect is 
by construction not included in standard thermal mod¬ 
els, generates a negative contribution to the second virial 
coefficient. As a result, for isospin-averaged observables 
where the isotensor channel has degeneracy (21 + 1) = 5, 
one finds an almost exact cancellation of the positive 
contribution due to /o(500) from the isotensor channel. 
This feature has been reported already in Ref. [13], and 
more recently in Ref. [14] , where similar conclusions were 
reached concerning the four-quark condensates and scalar 
susceptibilities of the hadron gas. 

As a net result of the cancellation, the combined effect 
of /o(500) and the isotensor-scalar channel on all isospin- 
averaged properties of the hadron gas is negligible. This 
is the basic result discussed in our study. 

A similar, albeit partial, cancellation occurs also in 
the pion-kaon S'-wave interaction. The attraction in the 
channel I = 1/2 is canceled by a repulsion occurring in 
the I = 3/2 channel. In the attractive I = 1/2 channel 
many studies find a kaon-like resonance, called a Kq (800) 
meson or n mm- 

The framework of the virial expansion in the context 
of hot hadronic gas was used in numerous works, includ¬ 
ing the nature of the repulsive terms mmm , meson 
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], thermal prediction from chiral perturba- 
3], or the treatment of the A(1232) reso- 


properties p 
tion theory 
nance [25H2T] . 

The paper is organized as follows: in Sec. [Tl]we briefly 
present the formalism and explicitly show the cancella¬ 
tion between the isoscalar and the isotensor channels in 


the S -wave (and also the D- wave) channels. In Sec. Ill 


we illustrate the cancellation with the trace of the energy- 
momentum tensor as a function of T, as well as for the 
abundance of pions. In Sec. |IV| we describe the simi¬ 
lar case of the kaonic resonance k. Section Ivl contains a 
discussion of two-particle correlation, where the cancella¬ 
tion mechanism does not occur and where /o(500) plays 
a relevant role [25]. Finally, in Sec. VI we present further 
discussion and conclusions. 


II. INTERACTING PION GAS 


For the case where the resonance is broad, or when 
we need to include interactions between stable particles 
which do not lead to resonances at all (such as repulsion), 
we should use the virial expansion. Then the 2 —>• 2 re¬ 
actions are incorporated according to the formalism of 
Dashen, Ma, Bernstein, and Rajaraman [29] 30J. The 
method, based on the virial expansion and valid for suf¬ 
ficiently low temperatures (T <C M, where M is the in¬ 
variant mass of the pair), uses as key ingredients the 
physical phase shifts which are well known for pion-pion 
scattering. Then 


d k (M) 


d5k(M) 

ndM 


(4) 


For simplicity of notation we focus on the case where 
the only degrees of freedom are pions. The partition 
function of the system up to second order in the virial 
expansion is given via the sum 


A. The formalism 


In modeling relativistic heavy-ion collisions, the study 
of hadron emission at freeze-out is successfully performed 
within thermal models (for reviews see, e.g., Refs. [8] !9j 
and references therein). In this framework, the free en¬ 
ergy of the hadron gas at temperature T, i.e., —ksT hi Z, 
is represented as the sum of contributions of all (stable 
and resonance) hadrons, whence 


In Z = hi Zf able + In Z( es , (1) 

k k 


where k indicates the particle species. In practice, one 
uses the list of existing particles from the PDG [I], 

In the limit where the decay widths of resonances are 
neglected, one has (for simplicity of notation we do not 
include chemical potentials) 
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where f k is the spin-isospin degeneracy factor, V is 
the volume, p is the momentum of the particle, the 
mass of the reson ance is denoted as M k , the energy is 
E p = sjp 2 + M|, and finally, the ± sign corresponds 
to fermions or bosons, respectively. Once In Z is deter¬ 
mined, all other thermodynamic quantities follow. 

As a better approximation for the partition function, 
one can take into account the finite widths of resonances. 
The quantity In Z k is replaced by the integral over the 
distribution function d k (M): 


In Zr = f k Vd k (M)dM 
Jo 


[ d \ In 

[i _ e-Ep/rl 

/ (27t) 3 



n -l 


( 3 ) 

For narrow resonances one can approximate d k (M) 
with a (non-relativistic or relativistic) normalized Breit- 
Wigner function peaked at M k . In the zero-width limit 
dfc(M) = 6(M — Mk) and Eq. ^ reduces to Eq. Q. 


In Z = In Z n +fu 
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where In Z n is the contribution from free pions evaluated 
according to Eq. ([2]), fjj = (2 1 + 1)(2 J + 1) is the spin- 
isospin degeneracy factor. In the case in which a series of 
narrow resonances labeled with k is present in the (/, J) 
channel, the derivative of the phase shift sharply peaks 
at the resonance positions according to 
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thus one recovers the contribution of the (/, J ) chan¬ 
nel to the formula of Eq. ([3]) in the Breit-Wigner limit. 
Moreover, for T'jj k —► 0 one obtains d8jj/{ndM) = 
Y^ k S(M — Mu t h), in agreement with Eq. ® [3D[. 

Quantum-mechanically, the quantity dS/JndM) has a 
very simple interpretation: it is the difference of the 
density of two-particle scattering states of the interact¬ 
ing and free systems. For completeness, we present the 
derivation of this well-known fact in the Appendix. 


B. The cancellation 

In Fig. [T] we report the experimental 7T7 t phase shifts as 
functions of the invariant mass up to 1 GeV and for dif¬ 
ferent values of (/, J), as parametrized in Refs. [6H7113T]. 
From panel (a) it is clearly visible that the channel (0, 0) 
is attractive (it is responsible for the emergence of the 
/o(500) pole), while the channel (2,0) is repulsive. The 
channel (1,1) is also attractive and corresponds to the 
prominent p meson. The isoscalar and isovector D -wave 
channels (J = 2) are reported in Fig. showing that 
the interaction strength in these channels is negligible. 

In Fig. [2] w e present the derivatives of the phase shifts 
from Fig. Ill multiplied with fu/n. We note from 
panel (a) that the resonance /o(500) does not lead to 
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FIG. 1. (color online) The parametrization of the pion-pion 
scattering phase shifts plotted vs invariant mass for various 
isospin-spin (JJ) channels, taken from Ref. [6l 171131 j. 


a pronounced peak, but rather to a smooth plateau with 
an integrable singularity close to the pion-pion thresh¬ 
old, which follows from kinematics. Importantly, the 
isotensor-scalar channel has a distribution, which up to 
M ~ 0.85 GeV is nearly a mirror reflection of the /o 
channel. Note that this is achieved with the multiplica¬ 
tion of the isotensor channel by the isospin degeneracy 
factor (2 1 + 1) =5, which occurs for isospin-averaged 
quantities (see Sec. Mfor other cases). As a result, the 
sum of the (0,0) and(2, 0) channels nearly vanishes until 
/o(980) takes over above M ~ 0.85 GeV. We note that 
the negative phase shifts in certain channels may be at¬ 
tributed to repulsion due to the finite size of hadrons, 
see, e.g. fI5 ]IT7H22] . 

Recall that chiral perturbation theory (yPT) predicts 
cancellation between the S -wave isoscalar and isotensor 
channels (for a historical review see, e.g., Ref. [ 22 ]). The 
5-wave scattering lengths a/o, related to the phase shifts 
near the threshold via Sjq = qajQ + o(q 2 ), with q = 
^/\pjA — read explicitly at lowest order %PT [33]: 


1 ni 7V 1 ^7T /—\ 

° 00 "32 nj*' a20 = ^7T (7) 


We note the opposite signs, hence cancellation, while 
the relative strength between the isoscalar and isoten¬ 
sor channels is —3.5. Inclusion of chiral corrections [24] 
brings the ratio to the value —4.95 ± 0.16, while a global 
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FIG. 2. (color online) The distributions 1/nfjj dSjj/dM , cor¬ 
responding to the phase shifts of Fig. [T] plotted as functions 
of the 7T7r invariant mass M. We note a nearly complete can¬ 
cellation for the S-wave interactions between the 1 = 0 and 
1 = 2 channels up to M ~ 0.85 GeV. 


analysis including dispersion relations made in Ref. |35] 
gives —4.79 ±0.55. These values are compatible with —5 
within the uncertainties. While yPT explains the can¬ 
cellation near the threshold, its occurrence until as far 
as M ~ 0.85 GeV is a quite remarkable fact which goes 
beyond standard yPT. 

In Fig. |2jb) we display the p-meson channel, which 
dominates the dynamics due to the fact that it is a pro¬ 
nounced resonance, and also carries a large spin-isospin 
degeneracy factor. In the remaining channels the inter¬ 
action is very small. In addition, in the D -wave channel 
there is also a partial cancellation between the (0,2) and 
(2, 2) channels at sufficiently low M. 


III. CANCELLATION OF / 0 (500) IN THERMAL 
MODELS 

As the cancellation between the attraction from 
/o(500) and the repulsion from the isotensor-scalar chan¬ 
nel occurs at the level of the distribution functions 
dk{M), it will generically persist in all isospin-averaged 
observables. The purpose of the examples in this section 
is to show, how much error one would make by solely 
including the /o(500) and neglecting the isotensor repul- 
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FIG. 3. (color online) Pion, p-meson, isoscalar-scalar, and 
isoscalar-tensor contributions to the volume density of the 
trace of the energy-momentum tensor divided by T 4 , plotted 
as functions of T. 

sion. 


A. Trace of the energy-momentum tensor 

In Fig. [3] we show the contribution to the trace of the 
energy-momentum tensor of the hadron gas from various 
channels: as expected, the contribution of a free pion 
gas is dominant at low T, while the contribution of the 
I = J = 1 channel (p meson) rises fast and overcomes 
that of free pions due to a large degeneracy factor. We 
note that the contribution from the attractive isoscalar- 
scalar channel (dashed line) is nearly canceled by the 
repulsive isotensor-scalar channel (dotted line). 

B. Abundance of pions 

The results of this subsection are of relevance to ther¬ 
mal modeling of particle production in ultra-relativistic 
heavy-ion collisions. In a simple version of this approach, 
hadrons (stable and resonances) are assumed to achieve 
thermal equilibrium, and later the resonances decay, feed¬ 
ing the observed yields of stable particles [32 38]. The 
basic outcome are thus the abundances of various hadron 
species, or their ratios, which allow to determine the val¬ 
ues of the thermal parameters at freeze-out. The con¬ 
tribution of light resonances is particularly important, 
therefore the treatment of / 0 ( 500) is relevant. 

In our study we use the SHARE nnum code to inves¬ 
tigate the effect of / o (500) on pion abundances. The 
purpose of this study is to see, how much pions are 
generated from the naive approach, where /o(500) is 
improperly included as a wide Breit-Wigner pole with 
M a = 484 MeV and T a / 2 = 255 MeV. In Fig. [2] we 
show the thermal model calculations carried out along 
the same freeze-out line as in Wl and for the tempera¬ 
ture of T = 156 MeV and all chemical potentials equal to 



FIG. 4. (color online) Relative contribution to the pion yield 
from /o(500), evaluated with SHARE [TPHUB] applied to A+A 
collisions at various collision energies. The thick lines corre¬ 
spond to the naive Breit-Wigner implementation of /o(500), 
while the thin lines indicate the correct inclusion of combined 
isoscalar-scalar and isoscalar-tensor channels. The points in¬ 
dicate the corresponding contributions at the LHC. 

zero at the LHC [571155] (the extra LHC points in Fig.[4j). 
We find that with the naive implementation, the relative 
feeding of the /o(500) to the pion abundances (i.e., the 
ratio of the number of pions originating from the /o(500) 
to all pions) would be at a level of up to 5%, which is of 
a noticeable size and in agreement with the conclusions 
of Ref. \iv , . 

However, in reality this contribution is not there due 
to the cancellation mechanism with the isotensor-scalar 
channel. The proper implementation of both the scalar- 
isoscalar and scalar-isotensor channels according to the 
phase-shift formula 0 for the density of states yields a 
very small negative contribution of the combined scalar- 
isoscalar and scalar-isotensor channels (at the level of 
—0.3%) to the pion yields, as indicated in Fig. [4] by the 
thin lines and the point below zero for the LHC. 

We note that the proper implementation described 
above leads to lower pion yields, thus higher model 
hadron to pion multiplicity ratios, at a level of a few 
percent. Of particular interest here are the kaon to pion 
ratio, /\ + / 7r+ , and the proton to pion ratio, p/n + . For 
the former, the cancellation mechanism results in higher 
values of the thermal model results. This may help to de¬ 
scribe the horn structure [55] at ,/snn = 7.6 GeV HD] , 
where models fall somewhat below the data HQ- The 
latter is related to the LHC proton-to-pion puzzle [4l1 . 
where the thermal model noticeably overpredicts p/n + . 
With the cancellation of the /o(500) contribution, the 
proton-to-pion puzzle becomes even stronger, opening 
more space for possible novel interpretations [42] . 

IV. THE CASE OF Ao(800) 

According to various works nuns], the attractive nK 
channel with I = 1/2 and J = 0 is capable of the gener¬ 
ation of a pole corresponding to the putative resonance 
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FIG. 5. (color online) The pion-kaon S-wave phase shifts (a) 
and their derivatives (b), multiplied with the isospin-spin de¬ 
generacy factors, plotted as functions of the invariant mass. 


A/) (800). This resonance has not been included in the 
summary table of the PDG, but it is naturally expected 
to exist as the isodoublet partner of the established res¬ 
onances /o(500) or / o (980). This wide resonance has a 
predicted mass of ~ 680 MeV and a very large decay 
width of ~ 550 MeV, which makes its final assessment 
quite difficult (see also the discussion in ‘Note on scalar 
mesons below 2 GeV’ in Ref. CD- 

For the pion-kaon system, we use the effective-range 
fits to the phase shifts given in Ref. dHl- The results 
are plotted in Fig. [5j where the phase shifts and their 
derivatives (multiplied with the isospin-spin degeneracy 
factors) for the channels (0,1/2) and (0,3/2) are pre¬ 
sented. Panel (b) also displays the sum (solid line), which 
exhibits the partial cancellation of the two contributions. 
Thus the conclusion here is similar, but less pronounced, 
compared to the case of the a meson. The I = 1/2 chan¬ 
nel is dominated by the vector channel J = 1, represented 
by the resonance A'* (892). 


V. CORRELATIONS 

Up to now we have considered only isospin-averaged 
observables, where the isospin degeneracy factor of 
(21 + 1) is present, leading to cancellation between the 


(0,0) and (2,0) channels. The cancellation does not oc¬ 
cur for correlation variables. Suppose, as an example, 
that we investigate the correlated 7r + 7r + pair production. 
Then, their source will be the isotensor channel only, and 
none will come from /o- 

An investigation of this type has been reported in 
Ref. [23, devoted to the analysis of the STAR Collab¬ 
oration data j44j on resonance production. In particular, 
for the 7T + 7r - pairs, the isospin Clebsch factors yield for 
the number of pairs 

n n+n -(M) = 3n p o(M) + n f o(M) + ^n (2 ,o)(M), (8) 

where the factor of 3 in front of the po contribution comes 
from the spin degeneracy. Since up to M ~ 0.85 GeV (cf. 
Fig. [2^,) the ratio n^ 2 ,o)(M)/nf 0 (M) ~ —1/5, the relative 
yield - ! of the isotensor channel compared to the isoscalar 
channel in Eq. Q is only about 10%, and there is no can¬ 
cellation. We thus clearly see the potential importance 
of the <r in studies of pion correlations [28] ■ 

VI. CONCLUSIONS 

In this work we have presented, in the framework of 
the virial expansion, that the contribution of the reso¬ 
nance /o(500) (alias a) to isospin-averaged observables 
(thermodynamic functions, pion yields) is nearly per¬ 
fectly canceled by the repulsion from the isotensor-scalar 
channel. The cancellation occurs from a “conspiracy” 
of the isospin degeneracy factor, and the physical values 
for the derivative of the phase shifts with respect to the 
invariant mass. Thus, the working strategy in thermal 
models of hadron production in relativistic heavy-ion col¬ 
lisions, which incorporate stable hadrons and resonances, 
is to leave out /o(500) from the sum over resonances. In¬ 
cluding it, would spuriously increase the pion yields at a 
level of a few percent. Alternatively, one may use phys¬ 
ical phase shifts to do more accurate studies, but then, 
of course, all relevant isospin channels must be incor¬ 
porated. The validity of the conclusion is for tempera¬ 
tures T smaller than the invariant mass of the hadronic 
pair. The lowest mass comes from the two-pion thresh¬ 
old, hence T < 300 MeV, which is a comfortable bound. 

A similar cancellation occurs for A/) (800) (alias k) with 
quantum numbers / = 1/2, J = 0, whose contribution 
is compensated by the I = 3/2, J = 0 channel. On the 
other hand, in correlation studies of pion pair production, 
there is no cancellation mechanism. 
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Appendix A: Derivation of the phase-shift formula 

For completeness, in this Appendix we quote the 
quantum-mechanical derivation of the phase-shift for¬ 
mula of Eq. ([3]). The relative radial wave function of 
a pair of scattered particles with angular momentum l, 
interacting with a central potential, has the asymptotic 
behavior ipi(r) ex sin[fcr — lir/2 + 5/], where k = |fc| is 
the length of the three-momentum, and Si is the phase 


shift. If we confine our system into a sphere of radius R, 
the condition kR — lir/2 + Si = rnr with n = 0,1, 2,... 
must be met, since ipi(r) has to vanish at the boundary. 
Analogously, in a free system kR — lir/2 = rif ree n. In the 
limit R —> 00 , upon subtraction, 

— = n - nfree■ (Al) 

IT 

Differentiation with respect to M yields immediately the 
interpretation, that the distribution d5i/(irdM) is equal 
to the difference of the density of states in M of the 
interacting and free systems un¬ 
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